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Introduction 

A recent  ARL  report  (Reference  1)  dealt  with  the  calculation  of  the 
mean  flow  of  a heated,  laminar,  water  boundary  layer  on  an  axisymmetric 
body.  In  the  development  of  the  pertinent  equations  a coordinate  trans- 
formation was  used  to  eliminate  explicit  dependence  on  the  coordinate  in 
the  flow  direction.  The  computed  velocity  and  temperature  profiles  were 
then  given  by  similarity  solutions.  The  detailed  laminar  boundry  layer 
characteristics  that  resulted  from  these  solutions  were  used  as  a means 
of  developing  a design  procedure  for  specifying  the  distribution  of  wall 
temperature  necessary  to  keep  the  flow  stable. 

Although  the  procedure  worked  well,  the  required  number  of 
parameters  made  the  interpolation  process  that  was  used  rather  cumbersome 
and  any  presentation  of  data  in  tabular  or  graphical  form  was  extensive. 
It  is  possible  with  a different  transformation  to  reduce  the  number  of 
parameters  by  one,  thus  reducing  the  complexity  of  the  computations  and 
ultimately  the  design  procedure. 

The  details  of  the  transformations  used  in  both  boundary  layer 
computation  methods  and  a description  of  the  interrelation  of  the  para- 
meters will  be  presented.  As  in  Reference  1,  the  effects  of  heating  on 
laminar  separation  and  transition  will  he  given. 

Development  of  Equations  and  Their  Solutions 

In  Reference  1,  for  steady  mean  flow  the  incompressible  boundary 


layer  equations  on  a body  of  revolution  were  presented  as: 
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3(p*r*u*)  3(p*r*v*)  _ 

3x*  3y* 


p*(u* 


v* 


3u*. 

3y*J 


dU  * . _ . 

5*V  “Si  + 5yi  <u‘  !^> 


(1) 


3y* 


3y*' 


3T*  3T*  ^ e*  3 3T* 

P*CP*<“*  Si  + v*  fi>  ' -p*“*V  -SI  + (k*  & 

The  boundary  conditions  for  these  equations  are: 

at  y*  = 0 u*  = v*  = 0 

as  v*  -*■  °°  u*  U * 

e 

T*  -*•  x * 


The  third  of  Equations  (1)  is  an  energy  balance  equation  in  which  buoyancy 
and  dissipation  by  friction  are  ignored.  The  stars  indicate  dimensional 
quantities.  In  dimensionless  terms  (see  List  of  Symbols),  Equations  (1) 
become: 


d(pru)  3 (prv)  _ 

3x  3y  “ v 


, 3u  3uv  _ ^Ue  ,13.  3u. 

p(u  3^  + V 3?}  ' pUe  + T W (U 


dU 


pc  (u  = -EPuU  J—  -r  - _ K—  -R- 

p 3x  3y  e dx  PRRe  °y  3y 


e + 1 


4-o.4?> 


Equations  (1)  and  (2)  do  not  include  transverse  curvature  terms; 


derivatives  of  r*(=  r*Q(x*)  + y*  cos  4>)  with  respect  to  y*  are 
small.  In  effect,  then,  r*  = rQ*(x*).  Also,  the  Eckert  number, 
evaluated  for  water  at  expected  conditions  in  the  free  stream  is 
quantity  and  the  term  in  which  it  appears  can  be  dropped. 


(2) 


i.e. , 

considered 
E,  when 
a small 


A stream  function  can  be  defined  as: 
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where  N Is  the  transformed  y coordinate. 


U * h 

N- 


<§-**>  <>y*- 


Using  Equations  (3)  and  (4)  along  with  the  definitions: 


. dU  * 
x*  e 

U * dx* 

e 


dr  * 

2 x*  o 

r * dx* 

o 


the  last  two  of  Equations  (1)  reduce  to: 


where  0 is: 


(pyF")  ' + M(^  - F’2)  + (M  t A +— *)  FF"  = 0 
PRC  ('W  + f'"  1-)  F0'  + (Pk0')’  = 0 


X*  - T * 

00 

X * - T * 

w 00 


The  boundary  conditions  become: 


at  N = 0 


F'  = F = 0 


0=1 


as  N -*•  <» 


F'  •*  1 


0 -*■  0 

Equations  (7)  were  solved  with  a slightly  altered  scheme  of 
Lowell  and  Reshotko  presented  in  Reference  2.  This  scheme  is  based 
on  the  numerical  integration  process  developed  by  Nachtscheim  and  Swigert 
(Reference  3).  The  procedure  that  was  adopted  was  to  assume  a series  of 
values  for  the  parameter  M,  A,  and  AT*  (=  Tw*-Tro*) , tabulate  the  calculated 
boundary  layer  characteristics,  and  use  these  results  in  a design  procedure. 
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If  Equations  (2)  are  transformed  according  to  the  Mangler-Levy-Lees 


transformation. 


d£  = U r dx 
e o 


dr)  = (Re/2^)^  U r dv, 


and  a reduced  stream  function. 


*(x,y)  = (2£/R  )**  f(5,n) 


and  the  parameter, 


,r  dU 
o _ 2£  e 

B " U dE, 
e 


are  introduced,  Equations  (2)  reduce  to: 


(PUf")'  +ff"  + 0(^  - f'2)  = 2£(f'f'  - f"f  ) 


f"  (Pkg')'  + Cfg'  = 2£Cp(f'g£  -g'f-). 

R 

g is  used  here  in  place  of  T.  If  the  functions  f and  g are 
assumed  to  be  independent  of  the  transformed  coordinate,  £,  the  following 
similar-type  homogeneous  equations  result: 

(PUf")'  + - f’2)  + ff"  = 0 


P_C  fg’  + (pkg')’  = 0. 

K p 

The  attendant  boundary  conditions  are: 
at  n - 0 


as  r) 


f = f'  = 0 


8 = T 


f’  - 1 


g 1 
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Although  Equations  (7)  and  (14)  were  derived  by  use  of  different  trans- 
formations, their  makeup  suggests  that  they  can  be  related  with  an  affine 
transformation  summarized  as  follows: 


Set  F = Af, 

n = bn. 

(15) 

Therefore, 

d _ d dn  d 

dN  dn  dN  B dn  ’ 

(16) 

F' 

= = AB  4^  = ABf ' , 

dN  dn 

(17) 

F"  = AB2f", 

(18) 

and 

F'  ” = AB3f"’  . 

(19) 

The  relation  between  g and  0 can  be  developed  with 


6’  - € - .(V  , dj. 


T * 


> ^ = b(- 


)g' 


and 


dN  T *-T  *'  dn  *-T  * 

w oo  w 00 


T * 

g"  = b"- 

w 00 


Making  these  substitutions  in  Equations  (7), 


+ A + 1.  . . 2„2 


(20) 


(21) 


(ABJ)  (pyf")  ' + (---'- y-'-x)  (AxP")ff"  + M(~  - A2B2f’2)  = 0 


prCp  (^^Hab)^ 


T * 

*)  fg'  *■  H2(^v.t  »>  (Pkg')'  = 0 (22) 

p w oo 


From  Equation  (17)  and  the  fact  that  in  both  systems  f'  = F'  = u/Ue> 


it  follows  that  AB  = 1.  Furthermore,  to  reduce  Equations  (22)  to 
Equations  (14), 
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A = 


B = 


M + X + 1 


M + X + 1 


I >S 


1 *5 


and 


, M_  - 2M 

M + X + 1 

D 


(23) 


(24) 


(25) 


Therefore,  the  M and  X parameters  can  be  combined  to  give  one 
parameter,  8.  X and  M are,  of  course  known  for  points  on  a given  body 
from  the  details  of  its  contour  and  the  pressure  distribution,  measured 
or  computed  by  way  of  the  Douglas-Neumann  procedure.  The  conclusion  is, 
of  course,  that  Equations  (14)  are  basically  the  equivalent  of  Equations 
(7)  and  can  be  used  to  generate  the  data  for  the  design  procedure  pre- 
viously mentioned.  In  fact,  solution  of  Equations  (14)  is  less  difficult. 

A finite  difference  procedure  due  to  H.  B.  Keller  [4]  was  used  in 
solving  Equations  (14).  The  procedure,  known  as  the  box  method,  was  deve- 
loped by  Keller  and  applied  successfully  by  Cebici  and  Smith  [5]  to  a wide 
variety  of  problems.  The  physical  properties  of  water,  C^,  p,  u and  k are 
a function  of  temperature  and  are  known  from  information  presented  in 
Reference  2. 


Results  of  Boundary  Layer  Computations 

The  boundary  layer  computations  were  performed  for  a parametric  study 
in  which  the  parameters  were  8 and  AT*  (“T  ^T^*) . In  all  cases  the 
free-stream  temperature,  T^*,  was  kept  at  60°  F.  A summary  of  the  results 
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of  this  parameter  study  is  found  in  Table  1.  Plots  of  these  data  are  also 
displayed  in  Figures  1-5. 

The  boundary  layer  thicknesses,  6*  and  0,  are  defined  in  the  sense 
of  two-dimensional  definitions, 

OO 

6*  = (l-p*u*)dy  (26) 

0 

CO 

0 = j p*u*(l-u*)dy*.  (27) 

0 

This  leads  eventually  to  slight  differences  when  comparing  with  thicknesses 
resulting  from  axisymmetric  considerations.  Little  difference  in  the  shape 
factor,  H (=  6*/0),  should  result,  however,  when  the  point  on  the  body  under 
consideration  has  a small  boundary  layer  thickness  compared  with  the  body 
radius.  For  considerations  of  a laminar  boundary  layer  this  is,  of  course, 
generally  true. 

The  quantities,  hg*  and  Hq,  given  in  Table  1 and  shown  in  Figures  1 
and  2,  must  be  transformed  to  get  6*  and  0 as  follows: 


as  critical  Reynolds  number  versus  temperature  difference.  The  A and 
AT  parameters  were  translated  into  the  shape  factor  H = 6*/0  so  that 

t 

in  each  case  R-*  . could  be  plotted  against  H.  The  curves  in 

o cnt 

Figure  6 are  the  result.  It  can  be  seen  that  the  correlation  is  quite  , 

good,  indicating  that  the  stability  of  the  laminar  boundary  layer  is 
strongly  dependent  upon  H,  regardless  of  whether  it  is  obtained  by 
favorable  pressure  gradient  or  by  heat.  A similar  correlation  is  reported 
in  Reference  7. 

This  correlation  provides  the  basis  for  the  development  of  the  type 
of  design  information  which  will  be  presented  here.  Although  based  on  two- 


dimensional  stability  information,  the  correlation  is  considered  valid  for 
the  axisymmetric  case  because  the  stability  equations,  under  the  assumption 
that  the  boundary  layer  thickness  is  small  compared  to  the  local  body  radius, 
are  the  same  for  both  cases. 


The  curve  in  Figure  6 for  the  variation  in  pressure  gradient  is  re- 


I 


k 


Determination  of  Body  Temperature  Distribution 

The  procedure  outlined  here  is  essentially  a simplified  version  of 
that  outlined  in  Reference  1.  Two  criteria  were  chosen  to  insure  the 
maintenance  of  laminar  flow:  (1)  the  provision  of  just  enough  heat  to 
keep  the  Reynolds  number  (based  on  displacement  thickness)  equal  to  the 
critical  Reynolds  number  and,  (2)  the  provision  of  enough  heat  so  that  the 
peak  critical  Reynolds  number  is  maintained.  These  are  referred  to  as 
"minimum  heat"  and  "maximum  heat"  conditions  respectively.  The  minimum 
heat  criterion  implies  that,  for  a particular  free-stream  velocity,  enough 
heat  is  added  to  make  the  operating  Reynolds  number  R^*,  equal  to  the 
critical  value.  This  should  insure  that  there  will  never  be  any  amplifi- 
cation of  waves  in  the  laminar  boundary  layer.  The  maximum  heat  criterion 
fixes  the  R,5*crit  at  its  maximum  value.  Whether  there  is  amplification 
depends  upon  the  free-stream  velocity  being  high  enough  to  have  the 

operating  Rx*  exceed  the  maximum  R-*  . . 

o 0 crit 

In  implementing  the  temperature  hunting  procedure,  the  data  in  Table  1 
can  be  filed  and  then  recovered  by  the  computer  for  interpolation  purposes. 
The  geometry  and  potential  flow  pressure  distribution  for  a body  will  be 
known  so  that  an  M and  a \ and,  consequently,  a 6 can  be  determined 
for  each  point  on  the  body  under  consideration.  This  will  be  designated  as 
Bq  for  a particular  body  point.  The  AT*  required  at  a body  station  can 
be  determined  by  applying  one  of  the  criteria  already  mentioned  and  inter- 
polating the  data  to  get  appropriate  quantities  corresponding  to  Bq. 
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In  the  minimum  heat  case  the  procedure  would  be  as  follows: 

(1)  Determine  an  Hq  versus  AT*  curve  corresponding  to  (3q  as 
illustrated  by  the  interpolation  procedure  pictured  in  Figure  3. 
The  curves  can  be  represented  by  spline  fits  and  values  are 
extracted  accordingly.  All  the  other  quantities  (r)^*,  qc>  T ) can 
be  similarly  handled  so  that  a curve  of  each  of  these  is  known  as 
a function  of  AT*  for  particular  8o's. 

(2)  First  choose  AT*  = 0.  The  n^*  corresponding  then  to 

AT*  = 0 and  can  be  used  to  determine  6*  for  a particular 

free-stream  velocity  using  Equation  (28). 

(3)  Calculate  R^*  by  means  of 

6 U * 


V 


e 

v * . 


(31) 


00 


This  is  the  operating  R^*. 

(4)  Enter  Figure  7 with  R^*  from  step  (3)  to  determine  a required 
H.  by  entering  with  the  operating  R^*  we  are  saying  that,  in  order 
for  this  to  be  the  ^,j*cr£t»  we  must  produce  a corresponding  value  of  H. 

(5)  The  H required  from  step  (4)  can  be  used  with  the  versus 

AT*  curve  of  step  (1)  to  determine  the  required  AT*.  This  required 
AT*  now  will  change  the  value  of  H^*,  originally  determined  in  step  (2) 
for  AT*  = 0. 

(6)  Repeat  ^tepr.  (2),  (3),  (4),  and  (5)  until  the  AT*  required 
converges  within  a desired  accuracy. 

(7)  Enter  the  q and  x versus  AT*  curves  established  for  8 

c o 

and  get  the  q and  T values  from  Equations  (30)  and  (31). 
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For  the  maximum  heat  case  the  somewhat  different  procedure  is  as 
follows : 

(1)  Assume  H = 2.29.  This  corresponds  approximately  to  the  point 

where  reaches  its  maximum.  This  is  not  a very  precise 

number  and  it  could  be  as  low  as  H = 2.2.  Adding  heat  beyond  the 

value  that  produces  the  maximum  Rr*  . will  result  in  R.* 

o crit  6 crit 

becoming  smaller  (see  Reference  4)  and  thus  be  counterproductive. 

(2)  Enter  the  curve  of  Hq  versus  AT*,  determined  in  step  (1)  of 

the  minimum  heat  procedure,  and  extract  the  AT*  required  for 

H = 2.29. 
o 

(3)  Extract  qc  and  for  the  AT*  of  step  (2)  from  curves  of 

qc  and  Tc  versus  AT*  also  determined  in  step  (1)  of  the  minimum 
heat  procedure  and  again  convert  to  q and  T via  Equations  (30) 
and  (31). 


Laminar  Separation 

An  attempt  was  made  to  obtain  limits  of  3 for  different  AT*  values 


beyond  which  laminar  separation  would  occur.  Using  the  criterion  that  the 


skin  friction  vanishes  at  the  point  of  separation,  curves  of  T versus 

c 

3 were  extrapolated  to  obtain  the  desired  limits.  The  computer  program  will 


not  calculate  a solution  to  the  boundary  layer  equations  at  T = 0 so  that 
extrapolation  is  necessary.  Figure  8 is  the  result  of  those  extrapolations. 


For  a given  local  temperature  difference,  laminar  separation  will  occur  for 
values  of  3 below  the  curve.  The  relatively  small  effect  of  temperature 
difference  is  apparent. 


r 
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As  a check  on  the  accuracy  of  the  calculated  separation  6,  an  attempt 

was  made  to  get  as  close  to  separation  as  possible  for  the  unheated  two- 

dimensional  case.  This  is  the  solution  for  Falkner-Skan  Flows  where  B 

is  now  the  Falkner-Skan  B.  These  results  were  then  compared  with  those 

presented  in  Reference  5 and  the  results  are  presented  in  Table  3.  The 

f"  values,  representing  the  slope  of  the  velocity  profile  at  the  wall, 

and  the  corresponding  B values  are  compared.  Richardson's  extrapolation 

discussed  in  Reference  5 was  used  to  get  the  ARL  values  of  f"  . 

w 

As  can  be  seen,  the  Smith  and  ARL  values  compare  favorably  near  laminar 
separation.  The  Keller-Cebeci  Value  of  B at  laminar  separation  does  not 
appear  to  plot  smoothly  with  the  rest  of  the  values  listed  and  thus  appears 
to  be  in  error.  A plot  of  f"  versus  B for  the  ARL  results  appears  in 
Figure  9.  Values  corresponding  to  reverse  flow  also  appear  in  this  plot. 


Transition 

The  prediction  of  transition  for  a given  temperature  distribution  on 

a body  follows  the  same  lines  as  presented  in  Reference  1.  The  basis  for 

the  method  used  here  is  the  plot  of  a hand  of  calculated  data  supplied 

through  the  courtesy  of  A.M.O.  Smith  (See  Figure  10).  The  band  marks  the 

a 

range  of  values  of  R versus  H that  were  obtained  from  e 

x trans 

stabilitity  calculations  performed  for  a variety  of  heated  and  unheated 
wedges.  It  would  be  logical  to  choose  the  lower  bound  of  this  band  as  the 
transition  criterion.  Values  for  such  a curve  are  given  in  tubular  form 
in  Table  4.  Translating  H into  8 for  this  transition  curve,  the  plots 
shown  in  Figure  11  result.  These  curves  indicate  the  transition  Reynolds 
number  that  would  correspond  to  particular  values  of  B and  AT*. 
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• As  stated  in  Reference  1,  the  curve  used  as  the  criterion  for  transition 

is  generally  optimistic,  that  is,  it  predicts  transition  at  a higher  Reynolds 
number  than  some  available  unheated  body  transition  data  would  indicate. 

This  is  evidenced  by  the  experimental  points  plotted  in  comparison  with  the 
transition  curve  in  Figure  12.  The  data  points  are  taken  from  Refererce  8. 
and  the  H values  for  these  data  were  obtained  from  the  information  in 
Table  1. 

In  work  not  reported  in  this  memorandum,  correlation  between  the 
Transition  Analysis  Program  System,  TAPS,  and  the  information  from  this 
similar  type  solution  indicates  that  the  H values  for  a given  temperature 
distribution  are  consistantly  higher  for  TAPS  by  a factor  of  about  1.03. 

This  correlation  is  with  respect  to  bodies  of  the  type  that  would  produce 
fairly  flat  pressure  distributions  and  which  would  profit  from  the  addition 
» of  heat  for  boundary  layer  stabilization.  A suggested  approach  for  esti- 

mating whether  or  not  transition  will  occur  is  to  replot  the  curve  of  Figure 
12  with  values  of  H reduced  by  1.03  and  then,  using  the  calculated  data 
of  Table  3 in  establishing  the  boundary  layer  characteristics,  one  would 
have  a curve  consistent  with  TAPS  and  experimental  transition  information. 

The  curves  of  Figure  11  could,  of  course,  be  similarly  altered. 
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Conclusions 

Techniques  developed  in  Reference  1 have  been  simplified  and  data  have 
been  generated  that  permit  one  to  estimate  in  a relatively  easy  manner  the 

temperature  distribution  necessary  to  stabilize  the  flow  over  an  axisymmetric 

» 

body.  Additionally,  these  techniques  and  data  provide  a means  for  deter- 
mining the  local  heat  flux,  skin  friction,  the  laminar  separation  point,  and 
the  point  of  transition  for  a heated  axisymmetric  body. 

All  these  estimates  can  be  made  by  using  the  values  of  the  various 
quantities  presented  in  the  tables  and  by  following  the  interpolative  pro- 
cedures and  criteria  that  have  been  outlined.  A simple  computer  program  to 
do  this  can  be  easily  written  or  one  can  interpolate  between  the  curves  that 


are  also  included  in  this  memorandum. 
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Computed  Laminar  Boundary  Layer  Data 
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Computed  Laminar  Boundary  Layer  Data 
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2.627405 
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60.0 
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0.512238 
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0. 103984 

87.657974 

-0 . loOCo 

60.0 

150.0 

1 .270015 

0.502926 

2. 525252 
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89.797714 

-0. 16000 

60.0 

150.0 

1 . 188935 

0.486326 

2.444728 
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Laminar  Boundary  Layer  Data 

T 

w 

V 

n0  H 

(°F)  (°F) 


qc 


— C.lcuQu  60.0 
-0.0 5000  60. 0 

.00000  60.0 
,0  500c  oC  • 0 
. lGCOO  60.0 
0.200 00  60.0 
0.30000  60.0 

0. 40000  60.0 

0.50000  60.0 


150.0 

160.0 
150.0 
150.0 
150.0 
150.0 
150.0 
150.0 
150.0 


1.027770 
0. 030634 
0.671094 
0.616824 
0.771875 
0.700824 
0.646428 
0.602951 
G. 567130 


0.446691 
0.420846 
0.399193 
0.380597 
0.364347 
0. 337094 
0.314944 
0.296456 
0.260713 


2.300852 
2.230622 
2. 182135 
2.146166 
2.116518 
2.079018 
2.052514 
2.033863 
2.020318 


0.299266 

0.352900 

0.399247 

0.440539 

0.478057 

0.544800 

0.603514 

0.656412 

0.704865 


100.450226 
104.606933 
107.893966 
1 10.618316 
112. 9464B  7 
116.784454 
119.880172 
122.473739 
124.704559 
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TABLE  2. 

Critical  Reynolds  Number  Variation  with  Shape  Factor 


6 . 

crit 

li 

1 

105.000 

3.80000 

2 

115.000 

3.5G000 

3 

130.000 

3.20000 

4 

150.000 

3.00000 

5 

180.000 

2.90000 

0 

232.000 

2.80000 

7 

326.000 

2.70000 

6 

485.000 

2.60000 

9 

675.000 

2.55000 

10 

1000.000 

2.51600 

1 1 

1500.000 

2.46750 

12 

2000. 0C0 

2.46600 

13 

3000. 0C0 

2.43750 

14 

4000.000 

2.41600 

15 

5000. 000 

2.40250 

lo 

6000.000 

2.36950 

1 7 

7000.000 

2. 37550 

18 

8000.000 

2. 36050 

19 

9000. 000 

2. 34400 

20 

10000.000 

2. 32  54  0 

21 

10500.000 

2. 31500 

22 

1 1000.000 

2.29750 

2 3 

11210. 000 

2.28  00  0 

tBIS 
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+.00064472, 

- . 0007  6458^ 

0 
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-.198838 
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-.195 
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-.195 

.0857052 
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-.190 
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-.180 
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-.180 
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.19078 

-.160 
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.239736 
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-.050 
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0 
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Transition  Reynolds  Number  Variation  with  Shape  Factor 
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trans 
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trans 
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09 
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09 
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09 
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08 
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08 
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08 
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07 
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07 
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07 
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07 

2.75C0C 

6.05486 
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07 
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06 
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06 
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06 
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5.33741 
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06 
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PRESSURE  GRADIENT  PARAMETER,  p 

Figure  3 - Shape  Factor  Variation  with  $ Parameter 
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Figure  8 - Laminar  Separation  Limit  Curve 
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Figure  9 - Variation  of  Pressure  Gradient  Parameter  with  Wall  Shear  Parameter  for  Falkner  - Skan  Flows 
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